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We present an analytical theory for the (K, if ')-protected core percolation problem. Each initially 
protected vertex transits to the susceptible state if it has less than K protected neighbors or if 
one of its susceptible neighbors has less than K' protected neighbors. The theory predicts (K, K')- 
protected core percolation to be a discontinuous transition for any type of uncorrelated networks and 
any K, K' > 2. We check these predictions by extensive computer simulations for the representative 
cases of K' = K. For random networks both the transition point and the {K, _R')-protected core size 
are predicted precisely by our theory. For low-dimensional lattices the transition point fluctuates 
considerably, but our theory still predicts precisely the {K, A')-protected core size once the transition 
occurs. The theory is also applied successfully to real-world complex networks. 

PACS numbers: 64.60.ah, 64.60.aq, 05.70.Fh 



Network percolation is underlying various phase tran- 
sition phenomena in physical and information systems 
[TJ [5]. The simplest percolation is the emergence of a 
giant connected component in the network (site or bond 
percolation) , which is deeply related to information prop- 
agation processes, such as disease spreading, opinion evo- 
lution, power-grid cascading failures, and resources trans- 
portation in biological systems [SHH]- The K-core of a 
network is the surviving subnetwork after all the ver- 
tices with less than K neighbors have been recursively 
removed from the network, if- core emergence is a per- 
colation transition [T0Hl2] . and it is strongly related to 
the phenomenon of particle jamming in glassy systems 
[Ol [Tl] and variable-state freezing in constraint satisfac- 
tion problems [T5] and Boolean dynamic networks (16j . 
Another network percolation process is the core forma- 
tion [TTUIQ) . The core of a network is obtained by a 
greedy leaf-removal procedure that repeatedly removing 
isolated vertices and vertices that are attached by a dan- 
gling edge. The existence of a percolating core causes 
great computational complexity in constraint optimiza- 
tion problems such as the minimal vertex cover problem 
PUI [21] and is also very significant for network control 
|22j . Network explosive percolations [53H1H]) are also in- 
tensively studied. 

Here we introduce and study a new network percola- 
tion problem, namely [K, if )-protected core percolation. 
In the model, besides the network connection pattern, 
each vertex also has a binary internal state (protected or 
susceptible) . Initially all the vertices are in the protected 
state. As time evolves, a protected vertex spontaneously 
transits to the susceptible state if it has less than K pro- 
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tected nearest neighbors. A protected vertex with K or 
more protected nearest neighbors is induced to the sus- 
ceptible state if at least one of its susceptible nearest 
neighbors has less than K' protected nearest neighbors. 
The [K, if ')-protected core is the subnetwork formed by 
all the surviving protected vertices and edges between 
them (Fig. [T] shows a simple example for K = K' = 2). 
For K' > 2, the (if, if ')-protected core percolation differs 
from the if-core percolation in that a susceptible vertex 
may induce the neighboring protected vertices into the 
susceptible state. 

(if, if ')-protected core percolation is qualitatively re- 
lated to many collective dynamical processes. One ex- 
ample is the frozen-core formation in Boolean satisfia- 
bility problems |291 I30j . where non- frozen vertices can 
induce its frozen neighbors into the non-frozen state un- 
der certain conditions (the whitening process [3TU33j ). 




FIG. 1: (color online). The (2, 2)-protected core (highlighted 
in red) of a small network. The (2, 2)-protected core is con- 
tained in the core (highlighted in green), which is again con- 
tained in the 2-core (highlighted in light blue). Protected and 
susceptible vertices are represented by white and black circles, 
respectively. 
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Another example is the glassy dynamics of kinetically 
constrained models, in which a spin in a certain state 
facilitates the flipping of its neighboring spins [34]. In 
inter-dependent networks, a collapsed vertex of one net- 
work causes the failure of the connected vertices in the 
other network [55H57] . sometimes resulting in a perco- 
lating damage cascading process. {K, iir')-protected core 
percolation also serves as a spreading model of a new 
product, with an adopted agent being able to persuade all 
his k non-adopted friends to adopt the product if k < K' . 
It may also be strongly related to jamming in transporta- 
tion systems [3H1 IMj • 

We present an analytical theory to describe {K, K')- 
protected core percolation. The theory predicts that 
the percolation transition is always discontinuous for any 
type of uncorrelated networks and any values of the con- 
trol parameters K,K' > 2. We confirm the theoretical 
predictions by performing simulations for the represen- 
tative cases of K' = K. 

Theory. — We consider a network formed by N vertices 
and M edges between pairs of vertices. If two vertices 
i and j are connected by an edge, then they are near- 
est neighbor to each other (self-connection and multiple 
edges between the same pair of vertices are not allowed) . 
The degree of a vertex is its number of nearest neighbors; 
and the degree profile P{k) quantifies the the fraction 
of vertices with each degree fc, with the property that 
X]fe>o ^(^) ~ 1- ^^'^ mean vertex degree of the network 
is c = 2M/N — J2k>o^-Pi^)- ^^'^ analytical theory 
the network is assumed to be uncorrelated. In such a ran- 
dom network the connection pattern is completely ran- 
dom, with only minimal degree-degree correlations be- 
tween nearest-neighboring vertices. We will study ran- 
dom regular (RR) networks, Erdos-Renyi (ER) networks 
and random scale- free (SF) networks [40l [41] . 

Each vertex of the network has two possible states, 
protected or susceptible. An initially protected vertex 
may spontaneously transit to the susceptible state or be 
induced to the susceptible state by its susceptible nearest 
neighbors. A non-empty (i^", i4r')-protected core of pro- 
tected vertices may form in the network after the state 



evolution reaches the final steady-state, each vertex of 
which having > K protected nearest neighbors and with 
each of the susceptible nearest neighbor having > K' 
protected nearest neighbors itself. For notational sim- 
plicity we discuss in the main text the representative 
cases of K' = K. The general cases of K' ^ K have 
the same qualitative properties [32]. In what follows we 
refer {K, i4r)-protected core simply as i^-protected core. 

Since the typical loop length of a random network 
scales logarithmically with vertex number TV, the net- 
work is locally tree-like [101 [33] ■ We therefore assume 
that given a vertex i being in the protected state, the 
states of all its nearest neighbors are mutually indepen- 
dent. This is a specific version of the Bethe-Peierls ap- 
proximation (also called the tree approximation), com- 
monly used in spin-glass theories, statistical inference, 
and complex network studies [3]-[51 [33] ■ Under this ap- 
proximation we then obtain an analytical expression for 



the normalized size rip-corc of the if-protected core: 
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In Eq. ([T]), a denotes the probability that a vertex j 
reached by following a randomly chosen edge (j,j) is in 
the susceptible state and having at most K — 2 protected 
nearest neighbors if not counting vertex i; and /3 is the 
probability that such a vertex j is in the susceptible state 
but having at least K — 1 protected nearest neighbors if 
not counting vertex i. Let us further define 7 as the prob- 
ability that such a vertex j is in the protected state and 
having exactly K protected nearest neighbors (including 
vertex i). The probability that vertex j has k nearest 
neighbors is denoted as Q{k). For an uncorrelated net- 
work Q{k) has the simple expression Q{k) = kP{k)/c 

Using the same Bethe-Peierls approximation, the fol- 
lowing self-consistent equations are obtained for the (cav- 
ity) probabilities a, /3 and 7: 
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The second term on the right hand side of Eq. ([2| is the 
probabihty that, when vertex j is still in the protected 
state, s > {K — 1) of its nearest neighbors (excluding 
vertex i) are also in this state. Even if vertex j in its 
protected state has > K protected nearest neighbors, it 
will still be induced into the susceptible state if it has 
a susceptible nearest neighbor which is surrounded by at 
most K—1 protected nearest neighbors. After j becomes 
susceptible, one or more of its protected nearest neigh- 
bors m (other than i) will spontaneously transit into the 
susceptible state if m originally has only K protected 
nearest neighbors. Then the number of protected near- 
est neighbors of j (not counting i) may be decreased to 
a value r < K—1. The last term of Eq. pi) is just the 
probability for this to happen. Equations ^ and Q are 
understood similarly. 

Equations (|3| and Q always have the trivial so- 
lution a — 1, p — J — 0, for which the normalized K- 
protected core size is np_corc = 0. This solution is al- 
ways locally stable ^35] , and it is the only solution if the 
network has too few edges. When edges are added to in- 
crease the mean vertex degree c, the degree profile P{k) 
will be changed. At certain critical value c — c* another 
stable solution suddenly appears. The value of a at c* 
is strictly less than 1, and np_corc jumps from to a 
value ?^p_coro that is strictly larger than [42 . The the- 
ory therefore predicts a discontinuous ii'-protected core 
percolation transition for any type of uncorrelated net- 
work and any value oi K > 2. The analytical theory 
also predicts that, for any K > 2 and independent of 
network type, fi-p-coro ~ "p-coro (c — c*)^/^ when c is 
only slightly beyond the critical value c* [H]. Such a 
scaling behavior was also observed in K-coie percolation 
and core percolation [TDHUl HH • 
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FIG. 2: (color online). Normalized iC-protected core size 
for random Erdos-Renyi networks. Solid lines are analytical 
predictions for K — 2 and K — 3 at N = oo, symbols are 
simulation results on three network samples with A'^ = 5 x 
10^, 10'*,2 X 10'^ vertices, respectively. 



ER networks. — The degree profile of the ER net- 
work follows the Poisson distribution when +oo, 
P(fc) = e~'^c'^ /k\ [33]. Figure [2] shows how np_core 
changes with mean vertex degree c. For K = 2 the the- 
ory predicts a percolation transition at c = c* w 3.9229, 
with a jump of np_core from to ".p_corc ~ 0.6181; for 
K = i the corresponding values are c* ~ 6.8426 and 
"p-coro ~ 0.7745. These predictions are fully confirmed 
by simulation results on single network instances. The 
2-core and core percolation transitions occur on ER net- 
works at c = 1 and c = e ~ 2.7183, respectively, and 
are both continuous [TH [15] . The 2-protected core per- 
colation transition is delayed to c w 3.9229 and becomes 
discontinuous but with critical scaling behavior. A dis- 
continuous transition was also discussed concerning dam- 
age propagation in inter-dependent networks |35l 136) . 

RR networks and lattices. — All the vertices in a RR 
network and a Z?-dimensional hypercubic lattice has the 
same degree fcoj a-nd the if-protected core contains the 
whole network when kg > K. If a randomly chosen frac- 
tion p of the edges are removed, the degree profile of 
the remaining network is P{k) ~ [ko\/kl{ko — k)l]{l — 
p)^ p^o-k ^ and the mean vertex degree is c = (1 — p)fco- 
The existence of a X-protected core then depends on c. 
For example, our theory predicts that n'^^^^^^ w 0.7725 
(for fco = 4) and « 0.7079 (for fco = 6) at the 

2-protected core percolation transition, with c* sa 3.0789 
and c* K, 3.3721, respectively. As demonstrated in Fig.[3j 
the theoretical predictions are in full agreement with sim- 
ulation results on single RR networks. 
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FIG. 3: (color online). Normalized 2-protected core size 
for RR networks with fco = 4 or fco = 6 (square symbols), 
the two-dimensional square lattice (cross symbols and dashed 
lines) and the three-dimensional cubic lattice (plus symbols 
and long-dashed lines). Each simulated network has A'^ = l(f 
vertices. The plus and cross symbols mark the transition 
points of 1000 independent instances of the square and cubic 
lattice system, respectively; the dashed and long-dashed lines 
are the results obtained on 20 independent instances of the 
lattice systems. The solid lines are analytical predictions. 
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FIG. 4: (color online). The average value of the 2-protected 
core percolation transition point c* for the D-dimensional hy- 
percubic lattice (square symbols). Each data point is an aver- 
age over 1600 network instances of A*' ~ 2 x lO'' vertices. Di- 
amonds are extrapolated simulation results to TV = -|-oo [42j . 
while circles are analytical predictions of c* for A'^ — +oo. 

The discontinuous 2-protected core percolation is also 
simulated on a Z3-dimensional hypercubic lattice with 
side length L and periodic boundary conditions {N — 
and /cq = 21?). For low dimensions the transition point c* 
is remarkably larger than the value predicted by the the- 
ory (see Figs. [3] and |4]), but the gap between the sample 
average of c* and the predicted value shrinks as D in- 
creases. The transition point c* fluctuates considerably 
for low dimensions (especially for D — 2,3) and also de- 
pends on N (for D < 7) [35j , and there is also no critical 
scaling behavior at c slightly beyond the transition point 
c* . Interestingly however, the value of ?^p_corc the 
percolation transition is located on the theoretical curve, 
and rip-corc follows the theoretical curve as c further in- 
creases (Fig. [3]). Much further work needs to be done 
on the if-protected core percolation in low-dimensional 
systems. 

Random SF networks. — Such networks are character- 
ized by a degree decay exponent A and a minimal degree 
fcniin. After a randomly chosen fraction p of the edges 
are removed, the degree profile of the remaining net- 
work is P{k) cx Em>max(/c„i„,fc) TO-^[m!/fc!(m - fc)!](l - 
p)^ p"^~^ , where the maximal degree fcmax is determined 
by NP{k^ax) ~ 1 [IS]- For single finite network instances 
the theoretical predictions on the X-protected core per- 
colation transition point and the normalized size np_coro 
(using only the P{k) as input) are in full agreement with 
simulation results [H] . As an example Fig. [5] shows the 
2-protected core percolation results for random SF net- 
works with A = 3 and /cmin — 4. There are strong finite- 
size effects: the transition point c* still increases with 
network size even for N > 10^. Such finite-size effects 
are due to the highly connected hub vertices, they be- 
come even more stronger as the exponent A decreases 



FIG. 5: (color online). Normalized 2-protected core size for 
random SF networks (A = 3 and fcmin = 4). Lines are the- 
oretical predictions with given vertex number A'^ and degree 
profile P{k), while symbols are simulation results on three 
independent instances with the same A*' and P{k). 
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Discussion. — This work demonstrated that the K- 
protected core percolation is a discontinuous transition 
for any type of uncorrelated networks and any K > 2. 
The same conclusion holds for the more general {K, K')- 
protected core percolation transition with K,K' > 2 
|42) . Our theory works both for infinite systems and sin- 
gle finite network instances. We have also applied the 
theory to 37 real-world networks, and the normalized 2- 
protected core sizes are correctly predicted for 32 of these 
networks [12] ■ 
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